arXi V: astro-ph/0112110 v2 11 Apr 2002 


submitted to The Astrophysical Journal 


Infrared Emission from Interstellar Dust. 

III. The Small Magellanic Cloud 

Aigen Li and B.T. Draine 

Princeton University Observatory, Peyton Hall, Princeton, NJ 08544, USA; 

agli@astro.princeton.edu, draine@astro.princeton.edu 

ABSTRACT 

The infrared (IR) emission from interstellar dust in the Small Magellanic Cloud 
(SMC) is modelled using a mixture of amorphous silicate and carbonaceous grains, 
including a population of polycyclic aromatic hydrocarbon (PAH) molecules. (1) It is 
shown that this dust model is able to reproduce the spectral energy distribution from 
near-IR to far-IR for the entire SMC Bar region, provided the PAH abundance in the 
SMC Bar region is very low. (2) The IR spectrum of the SMC Bl^l molecular cloud 
can also be reproduced by our dust model provided the PAH abundance is increased 
relative to the overall SMC Bar. The PAHs in SMC Bl^l incorporate ~3% of the SMC 
C abundance, compared to < 0.4% in the SMC Bar. (3) The spectrum of SMC Bl#l 
is best reproduced if the PAH mixture has intrinsic IR band strengths which differ 
from the band strengths which best fit Milky Way PAH mixtures. The variation in the 
PAH IR band strengths would imply different PAH mixtures, presumably a consequence 
of differing metallicity or environmental conditions. Other possibilities such as super¬ 
hydrogenation of PAHs and softening of the starlight spectrum are also discussed. 

Subject headings: dust, extinction — infrared: ISM: lines and bands — galaxies: indi¬ 
vidual: SMC 


1. Introduction 

Due to its low metallicity (~1/10 of that in the Milky Way [Dufour 1984]) and low dust-to-gas 
ratio (over 10 times lower than in the Milky Way [Bouchet et al. 1985]), the Small Magellanic 
Cloud (SMC) is often considered as a prototype for the interstellar matter in high-redshift galaxies 
at early stages of chemical enrichment. The analogy is strengthened by the similarities of its typical 
extinction curve - a nearly linear rise with (A is wavelength) and the absence of the 2175A hump 
(Lequeux et al. 1982; Prevot et al. 1984; Bouchet et al. 1985; Thompson et al. 1988; Rodrigues 
et al. 1997; Gordon &: Clayton 1998) - to wavelength-dependent extinction inferred for starburst 
galaxies (Gordon, Calzetti, & Witt 1997) and damped Lya systems (Pei, Fall, & Bechtold 1991). 
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Infrared (IR) emission from dust in the SMC has been seen by the Infrared Astronomical 
Satellite (IRAS) (Schwering & Israel 1990; Sauvage, Thuan, & Vigroux 1990), the Diffuse Infrared 
Background Experiment (DIRBE) instrument on the Cosmic Background Explorer (COBE) satellite 
(Stanimirovic et al. 2000), and the Infrared Space Observatory (ISO) (Reach et al. 2000). In the 
present work we use these observations to test the grain model proposed previously to account 
for the observed extinction curve for dust in the general SMC Bar region (Weingartner & Draine 
2001a). 

Of particular interest is the detection by Reach et al. (2000) of the 6.2, 7.7, 8.6, 11.3, and 
12.7/im emission bands toward the quiescent molecular cloud SMC Bl^l. These bands, generally 
attributed to polycyclic aromatic hydrocarbons (PAHs) (Leger &: Puget 1984; Allamandola, Tielens, 
& Barker 1985), had not previously been observed in the SMC. 

The molecular cloud SMC Bl#l ([a,djigso = [0'^43™42®4, —73°35'10"]) has a 11.3//m/7.7;um 
emission band ratio which differs significantly from that for Milky Way regions, such as the p Oph 
reflection nebula (which one might expect SMC to resemble [see Reach et al. 2000]); the 

observed band ratio also falls outside the range which appears to be “allowed” by the PAH model 
of Li & Draine 2001 (see Eigures 16,17 of Draine & Li 2001). The observed spectrum for SMC Bl^^^l 
therefore presents an important test which may reveal shortcomings of the current PAH emission 
models. The unusual band ratio might be indicative of the physical conditions in SMC or 

it might point to systematic differences between the PAH mixtures in the Milky Way and SMC. 

In this work we extend the interstellar dust emission model developed for the diffuse interstellar 
medium (ISM) of the Milky Way (Li & Draine 2001; hereafter LDOl) to the SMC, with emphasis 
on the mid-IR spectrum of SMC Bl^l. We seek (1) to test the PAH IR emission model (Draine 
& Li 2001; LDOl); (2) to infer the SMC PAH properties such as size distributions, charging, 
hydrogenation, and intrinsic IR band strengths; (3) to infer the physical conditions in the SMC 
and the SMC Bl^l cloud and the dust properties in these environments. Our approach is outlined 
in §2.1 and §2.2. We first model the spectral energy distribution for the SMC as a whole in §3. 
Model results for the SMC Bl^l molecular cloud are presented in §4.2-§4.5. In §4.6 we discuss 
the 2175A extinction hump predicted from the present PAH model. Our principal conclusions are 
summarized in §5. 


2. Modelling the Dust IR Emission 
2.1. Dust Model 

Eollowing LDOl, we model the SMC Bar and the SMC Bl^l molecular cloud IR emission 
as due to a mixture of amorphous silicate grains and carbonaceous grains. We assume that the 
carbonaceous grain population extends from grains with graphitic properties at radii a > 0.01;um, 
down to particles with PAH-like properties at very small sizes. 



3 


Dust grains in the SMC and SMC Bl^^l are assumed to be heated by a radiation field equal 
to the Mathis, Mezger, & Panagia (1983; hereafter MMP) estimate of the solar neighbourhood 
interstellar radiation field (hereafter ISRF),^ multiplied by a factor U. The IR emissivity per unit 
solid angle per H nucleon from a mixture of silicate and carbonaceous grains (including both neutral 
and ionized PAHs) is 
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where nn is the hydrogen number density; dnca,rh/da and dusn/da are the size distribution func¬ 
tions for carbonaceous and silicate grains, respectively; C'abs'^^’ ^abs absorp¬ 

tion cross sections for neutral and ionized carbonaceous grains (including PAHs) and silicate 
grains, respectively (see LDOl); Bx{T) is the Planck function; dP^^^'°° {U,a,T), dP^^^'°^ {U,a,T), 
and dP^'^^{U,a,T) are, respectively, the probabilities that the vibrational temperature will be in 
[T, T + dT] for neutral and charged PAHs and silicate grains illuminated by starlight intensity U 
{dP/dT approaches a 6 function for large grains); and (f)ion{U, a) is the probability of finding a PAH 
molecule of radius a in a non-zero charge state. 


For a given radiation intensity U, we calculate the temperature distribution functions dP/dT 
employing the “thermal-discrete” method which treats both heating (photon absorption) and cool¬ 
ing (photon emission) as discrete transitions, using a thermal approximation for the downward 
transition probabilities (Draine & Li 2001). We have calculated P(U,a,T) for neutral and ionized 
PAHs and silicate grains for 0.1 < U < 10^ at intervals of 0.25 in logj^QH. 

Except for single optically-thin clouds, the region observed will generally include material 
exposed to a range of radiation intensities. We follow Dale et al. (2001) and assume the radiation 
intensity U to have a power-law distribution of the general form 

^ _ (l-,3)ivr , ^ 

ITT ~ Trl—fS jjl—13'^ ’ ^ ^ Umax ; H ^ i \^) 

Pmax — c/jnin 
^tot 

= , for H < [/n,ax , /? = 1 ; (3) 

ln(L'niax/^min) 

where C/min and Hmax are, respectively, the lower and upper cutoff of the radiation intensity U, and 
is the total hydrogen column density. Lacking knowledge of the detailed geometry of the dust 


^The SMC Bar starlight spectrum may be harder than the MMP ISRF since the SMC has a higher surface density 
of O and B stars (Lequeux 1989). On the other hand, grains within the SMC Bl^l cloud are exposed to a softer 
radiation held due to reddening by dust in the cloud (see §4.4). However, in comparison with other uncertainties, 
these should be considered as second-order effects. We will see in §4.4 that reddening of starlight has only a minor 
effect on the PAH mid-IR emission band ratios. 
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and illuminating stars, we will assume the different illumination levels to be randomly distributed 
along the line-of-sight. The emergent IR intensity is then 


1 - exp [-A^^°tSabs(A) 
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where Sabs (A) is the total absorption^ cross section per H nucleon for the dust model. The integrated 
IR power per area is N^^Po{U), where po is the absorbed power per H nucleon for dust exposed to 
starlight with U = 1, and the mean starlight intensity (U) is 
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2.2. Model Parameters 

We have three sets of parameters to be specihed or constrained: 

1. The column density A^h of gas averaged over the beam. 

2. The parameters {he, ao, and a) determining the abundance and size distribution of the “log¬ 
normal” population which contributes most of the ultrasmall carbonaceous particles (PAHs). 

3. Environmental properties including the starlight intensity U (determined by l/min; Uma.x, and 
(3), the electron density Ug, and the gas temperature Tgas. The heating and cooling of dust 
grains are determined by U. The PAH ionization fraction (pion depends on Vjn^,, and Tgas- 

For dust in the general SMC Bar we will adopt the grain size distribution obtained by Wein- 
gartner & Draine (2001a) by fitting the extinction curve toward the star AvZ 398 in the SMC Bar 
(Gordon & Clayton 1998). We will see that the PAHs in the SMC Bar make a negligible contribu¬ 
tion to the IR emission. The SMC Bar models therefore depend on just 4 adjustable parameters: 
fAmin, Cmax, /?, and 

As noted by Reach et al. (2000), SMC Bl^l may be considered to be an analog of the p Oph 
molecular cloud. The dust in the p Oph cloud has Ry 4.2, so we will assume that the dust 
grains in SMC Bl^l have the size distribution estimated by Weingartner & Draine (2001a) for 
Milky Way dust with Ry = 4.0 (6c = 20ppm, Case A), except (1) with abundances relative to H 


the radiation field in the emitting region is approximately isotropic, scattering does not alter the intensity. 
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reduced by a factor 10, the SMC gas-to-dust ratio relative to the local Milky Way value (Bouchet 
et al. 1985), and (2) we allow ourselves the freedom to adjust the parameters be, ao, and a. 

We adopt a single log-normal size distribution for the PAHs in the SMC Bl^l cloud, charac¬ 
terized by three parameters: oq, ct, and be; uo and a respectively determine the peak location and 
the width of the log-normal distribution, and be is the total amount of C atoms relative to H locked 
up in PAHs. Small PAHs (with < 100 carbon atoms) are expected to be planar (see Appendix A 
in Draine & Li 2001). The term “PAH radius” used in this and related papers refers to the radius 
a of a spherical grain with the same carbon density as graphite (2.24gcm“^) and containing the 
same number of carbon atoms A^e'- a = 1.288N^‘^^ A. 

We adopt nu ~ 10^cm“^, and Tgas = 15K for the SMC Bl^^^l molecular cloud as estimated 
by Lequeux et al. (1994) based on observed CO line ratios. After allowance for C in the dust 
grains, the gas phase C abundance is [C/H]gas ~ 3 x 10“®. In diffuse and translucent regions with 
^ 3, we expect C and other gas-phase metals (Mg, Si, S, Fe) to be photoionized, contributing a 
fractional ionization ~ 4 x 10“^ in SMC gas. We consider two ionization models: low ng and high 
Hg. For the low-ng model we take rig = 0.1 cm“^ {nglnn ~ 1 x 10“^). For the high-Ue model we 
take Ug = 1.0cm“^ (ne/nn ~ 1 x 10“^). For the SMC Bar, we take Ug = 0.1 cm“^. But this is not 
critical since the amount of PAHs in the SMC Bar is rather small and the PAH mid-IR emission is 
negligibly low comparing to the far-IR emission (see §3). 

For a given Ug, our model for SMC Bl^^tl is then determined by 7 parameters: t/min; L^max; P; 
the dust column cx the fraction of the dust absorption due to PAHs oc be; and the parameters 
ao, CT characterizing the PAH size distribution. We will assume a = 0.4 for the width of the PAH 
size distribution as determined for the diffuse ISM PAHs (LDOI). Thus we are left with Apar = 6 
adjustable parameters: 17min, T/max, P, be, and oq. 


3. IR Emission from the SMC Bar and Eastern Wing 
3.1. Observational Constraints 

We first consider the emission from the SMC as a whole. In Figure 1 we plot the surface 
brightness measured by DIRBE for a 6.25 deg^ region containing the optical bar and the Eastern 
Wing (Stanimirovic et al. 2000). The overall spectrum, unlike that of SMC BI^t^I (see §4), peaks 
at A ~ 100//m with a local minimum at A ~ 12/im. 

Stanimirovic et al. (1999) have mapped the HI 21cm emission from the SMC. The total HI 
mass, corrected for self-absorption, is ~ 4.2 x 10 ®Mq (Stanimirovic et al. 2000). The 6.25 deg^ 
region considered here contains ~ 2.9 x 10 ®Mq of HI (S. Stanimirovic 2001, private communication), 
corresponding to a mean HI column density (AHi)obs ~ 5.5 x 10^^ cm“^, consistent with (Ahi) ~ 
5.0 X 10^^ cm“^ determined by Lyman a observations of 21 sightlines (Bouchet et al. 1985). Based 
on a FUSE {Far Ultraviolet Spectroscopic Explorer) survey of H 2 along 26 lines of sight, Tumlinson 
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et al. (2002) found the H 2 /HI mass ratio M(H 2 )/M(HI) < 0.5%. However, Israel (1997) suggest a 
global ratio M(H 2 )/M(HI) w 20%. We take (AiH)obs ± ctr = (5.5 ± 0.5) x 10^^ cm“^ for the region 
considered here. 

For a model with Aipar adjustable parameters, the goodness of fit is measured by 

^ ([A^Hlmod - (A^H)obs)Vo'H + ELi ([^-^A]mod - [A/A]obs)^ ! . . 

d.O.f. .^par + 7 

measuring the difference between model and observations for the 6.25 deg^ field studied by Sta- 
nimirovic et al. (2000). We use observed 12, 25, 60, 100, 140, and 240/im DIRBE band surface 
brightnesses [/A]obs and uncertainties ai from Stanimirovic et al. (2000). 


3.2. Model Results 

Grains over such a wide area will be illuminated by a range of radiation fields. We model the 
SMC DIRBE spectrum assuming a power law distribution for the illuminating radiation (see eq. 
[2]), using the HI column density (iVH)obs as an additional constraint. As discussed in §2.2, we are 
left with Npar = 4 adjustable parameters: I/min, t^max, [3, and 

Our best-fit model (see Table 1) provides a good fit to the entire spectrum with emission from 
a mixture of silicate and carbonaceous grains with ~ 5.4 x 10^^ cm“^, illuminated by starlight 
with t/min = 10“^'°, [/max = 10^'^®, (3 = 1.8 (see Figure la). The small value of x^/d.o.f. = 0.2 
suggests that some of the observational uncertainties may have been overestimated. 

Our best model has a dust-weighted mean radiation intensity ([/) ~ 1.9. The observed mean 
surface brightness in the 1.25^m DIRBE band averaged over the 6.25 deg^ held corresponds to 
only 0.76 of the intensity in the Mathis, Mezger, & Panagia (1983) radiation held, but we expect 
([/) to be larger since the dust is preferentially located in regions of high stellar density.^ Other 
distributions of radiation intensity are also able to provide a fairly good ht to the observational 
data, but only if [/min ^ 1 and [/max ^ 50 (see Table 1). 

We have also modelled the SMC data using a single [/, i.e., assuming that the dust is heated 
by a uniform radiation held with intensity U. The best-htting model is given by [/ = 10^'^ and 
w 2.8 X 10^®cm“^. Although this model provides a fairly good ht to the overall spectrum, 
the derived hydrogen column density is smaller than measured by a factor of ~ 20, giving 
X^/d.o.f. = 10.2. This model is rejected: the single U model has overestimated the mean radiation 
held and correspondingly underestimated 


®An enhanced radiation field 10 — 100 times the Galactic mean value) is indicated by the observed rotational 
excitation of H 2 in the SMC diffuse ISM (Tumlinson et al. 2002). This may be due to the fact that the massive stars 
and H 2 which FUSE observed are located near regions of higher-than-average starlight. 
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Table 1; Models for SMC Bar and SMC Bl^l. 


Region 

note 

ne^ 

(cm-3) 

PAH band 
strength^ 

ao 

(A) 

a 

6C 

(ppm) 

^min 

(MMP) 

t^max 

(MMP) 

0 

{U) 

(MMP) 

iVg’i 

(10^1 cm-2) 

x“/ 

d.o.f. 

(mag) 

■^^2175 A 
(mag) 

Bar 

best 

0.01 

MW 

- 

- 

0 

IQ-l U 


1.8 

1.9 

5.4 

0.2 

0.33 

- 

Bar 

Cmi„ T 

0.01 

” 

- 

- 

” 

10-0.5 

102.75 

2 

2.4 

4.8 

0.9 

0.30 

- 

Bar 

t^max i- 

0.01 

” 

- 

- 

” 

10-10 

102.0 

1.6 

2.3 

4.8 

1.0 

0.30 

- 

Bar 

/3 i 

0.01 

” 

- 

- 

” 

10-10 

IOI .75 

1.5 

2.4 

4.9 

1.4 

0.30 

- 

Bar 

0 T 

0.01 

” 

- 

- 

” 

IOO .25 

103.0 

2.5 

5.1 

2.4 

5.6 

0.15 

- 

Bar 

single U 

0.01 

” 

- 

- 

” 

IOI .5 

IOI .5 

- 

32 

0.28 

10.2 

0.02 

- 

Bl#l 

best 

1.0 

SMC 

2.8 

0.4 

1.5 

10-0.5 

102.5 

2 

2.2 

42 

1.7 

2.05 

0.5 

Bl#l 


0.1 

MW 

3 

0.4 

1.5 

10-0.5 

102.5 

2 

2.2 

42 

5.3 

2.05 

0.5 

Bl#l 


1.0 

MW 

” 

” 

” 

” 

” 

” 

” 

” 

3.1 

” 


Bl#l 


0.1 

Bq.(12) 



” 

” 


” 

” 

” 

3.3 

” 


Bl#l 


1.0 

Bq.(12) 



” 

” 


” 

” 

” 

1.9 

” 

” 

Bl#l 


0.1 

SMC 



” 

” 


” 

” 

” 

2.4 

” 

” 

Bl#l 


1.0 

SMC 



” 

” 


” 

” 

” 

1.8 

” 

” 

Bl#l 


0.1 

SMC 

2.8 

0.4 

” 

” 

” 

” 

” 

” 

2.1 

” 

” 


“Electron density Ue used to calculate the PAH ionization fraction <f>io-n, assuming Tgas = 15K. 

^'“MW” refers to the PAH band strengths adopted for the Milky Way diffuse ISM (LDOl) for which the 6.2, 7.7, and 
8.6/rm bands have been enhanced (relative to laboratory values) by factors of 3, 2, and 2, respectively; and the H/C 
ratio is taken as Eq.(ll); “Eq.(12)” refers to the same band strengths as those of LDOl but with a higher H/C (see 
Eq.[12]; “SMC” refers to the laboratory measured band strengths except the 6.2/im band is enhanced by a factor 1.5 
(with Eq.[ll]). 

“We adopt Ay = 4.9 x 10“^®, 6.2 x 10“^® magcm^ respectively for the SMC Bl#l molecular cloud and the SMC 
Bar (Bouchet et al. 1985; Martin, Maurice, & Lequeux 1989). 

3.3. The PAH Abundance in the SMC Bar 

The SMC Bar dust model described above consumes [C/H]dust ~ 15 ppm, and [Si/H]dust ~ 
12 ppm. Our SMC Bar grain model has a very low abundance of PAHs; the a < IsA carbonaceous 
grains comprise only [C/H]pah ~ 0.2 ppm, consistent with nondetection of the 2175A extinction 
hump (Weingartner & Draine 2001a). 

An upper limit on the PAH abundance can be obtained by comparing the DIRBE 12;um and 
25/tm photometry with theoretical IR spectra calculated for PAH molecules illuminated by the 
above derived radiation field. We consider PAHs of a single size a. For each size a we obtain an 
upper limit on [C/H]pah in PAHs such that this abundance of PAHs would not - by itself ~ exceed 
the DIRBE 12^m or 25/um photometry (i.e., we ignore possible additional emission from silicate 
grains and larger carbonaceous grains). In Figure 2 we show the resulting upper limit [C/H]pah as 
a function of PAH size a. It is seen that the DIRBE 25//m band places a stronger constraint than 
the 12^m band for PAHs > IDA; whereas for smaller PAHs (< IDA) the 12/im photometry provides 
a stronger constraint. Figure 2 also indicates that a ~ 12A grains are the best 12^m emitter: the 
upper limits [C/H]pah increase as grains either increase or decrease their sizes. In summary, we 
conclude that the DIRBE 12/im and 25/im photometry places an upper limit of [C/H]pah ^ 2 ppm 
on the abundances of PAHs (a < 15A). This upper bound is much less stringent than the upper 
limit [C/H]pah ^ 0.2 ppm implied by the absence of the 2175A extinction hump toward AvZ398. 

The paucity of PAH molecules in the SMC Bar, together with the fact that the average 
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Fig. 1.— Comparison of the model (solid line) to the observed emission from the SMC obtained by COBE/DIRBE 
(diamonds) and IRAS (squares) averaged over a 6.25 deg^ region including the optical bar and the Eastern Wing 


(Stanimirovic et al. 2000). Triangles show the model spectrum convolved with the DIRBE filters. Stellar radiation 


(dot-dashed line) dominates for A < G/rm. Grains are illuminated by a range of radiation intensities dNu/dU oc U 
10"^° <U < 10^'^®, with « 5.4 X 10^^ cm-^ 


extinction curve of the SMC has no observable 2175A hump, is consistent with the hypothesis that 
PAHs are the carriers of the 2175A feature. SIRTF 5 — 25/im observations of SMC regions where 
the 2175A feature is either very weak or absent could strongly test our dust model, which predicts 
minimal 7.7 and ll.S/tm emission features from such regions (as in Figure 1). 


4. IR Emission from the SMC Bl^^l Cloud 

Having seen that our dust model is able to closely reproduce the emission from the SMC as a 
whole, we now consider the emission from the SMC Bl^l quiescent molecular cloud. 


4.1. Observational Constraints 

The goodness of fit for SMC Bl^l is measured by 

_ El=iWi{[XhUd-[XhUs? /af 

d.o.f. Mp,, + IF, ^ ^ 

where A^par is the number of adjustable parameters, [A/Ajmod is the model spectrum (see Eqs.[l,4]), 
and [A/yJobs is the observed spectrum. 

We use the 5-15^m ISOCAM spectrum (Reach et al. 2000). For the 150 wavelengths in the 
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Nc 
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Fig. 2.— Upper limits (solid line) on the abundances of PAHs in the SMC Bar placed by the DIRBE 12pm 
and 25pm band photometry and the 2175A hump extinction (dot-dashed line). Overall, the DIRBE photometry 
limits the amount of PAHs in the SMC Bar to be [C/H]pah ^ 2 ppm (heavy solid line). The absence of the 2175A 
extinction hump provides the strongest constraint: [C/H]pah iS 0.2 ppm (Weingartner & Draine 2001a). The upper 
axis indicates the number of C atoms Nc in a PAH molecule. 


ISO spectrum, the weights Wi are 

taken to be 




/ 

1, 

for |A - Aa| < 7 q:A„ 

(a = 1,2,... 

5) 

Wi = < 

1/3, 

for |A - Aa| > 7 q,A„ 

(a = 1,2,... 

5) 


1/6, 

for A < 6.02//m or 

A > 13.38/tm; 



( 10 ) 


where Aq and 7 q,Aq are the peak wavelength and FWHM of the a-th PAH feature (see Table 2). 
and (Tj is the l-cr uncertainty in [A/Ajobs at wavelength Ap 

In addition, we use the IRAS 12, 25, 60, and 100/im intensities, with cjj = 0.1(A/A)obs (Schw- 
ering &: Israel 1990); and the DIRBE 140 and 240;um photometry (Hauser et al. 1998; Schlegel, 
Finkbeiner, & Davis 1998) after an empirical factor of 1.6 correction for beam dilution^, with 
cjj = 0.2(A//a)p We assign Wi = 1 for the 4 IRAS and 2 DIRBE bands, giving ^ lFj=i = 93. 


^The DIRBE beam size (0.7° x 0.7°) is much larger than the angular size of the SMC Bl^l cloud (~ 50” x 40”, 
physical size ~ 10pc), but there are a number of other emitting clouds in the DIRBE beam. To empirically correct 
for beam dilution, we multiply the DIRBE 60, 100, 140, and 240pm surface brightnesses by a factor of 1.6, bringing 
the DIRBE 60 and 100pm data into agreement with IRAS (see Figure 4). This provides our best estimate for the 
emission at 140 and 240pm. 
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4.2. Models with Milky-Way PAH Band Strengths 

For our initial model, we use IR cross sections C'abs(a, A) for neutral and charged PAHs following 
LDOl. We use the band strengths f Sa(A)dA~^ recommended by LDOl for the Milky Way (MW) 
PAH mixture where <S'q-(A) is the absorption cross section per C atom, but we adjust central 
wavelengths Aq, and FWHMs 7 q,Aq, of the mid-IR features to match the observed spectrum. Table 
2 shows the canonical MW band positions and FWHM (from LDOl) and the values appropriate for 
SMC Bl^l (see also Table 1 in Reach et al. 2000). Most of the changes are small, but the 11.3 and 
12.7/im features in SMC Bl#l are approximately twice as wide as for MW PAH mixtures, while 
the 7.7/im feature is ~25% narrower than in the MW. 


Table 2: Prude Profile Parameters for PAHs 


a 

MW'^ SMC^ 

7aAQ(^m)'^ 

MW^ SMC® 

MW'^ 

C 

y 

SMCf 

1 

6.20 

6.26 

0.20 

0.24 

3 

1.5 

2 

7.70 

7.65 

0.70 

0.53 

2 

1 

3 

8.60 

8.48 

0.40 

0.36 

2 

1 

4 

11.30 

11.34 

0.20 

0.54 

1 

1 

5 

12.70 

12.80 

0.30 

0.58 

1 

1 


“Central wavelength. 

*'FWHM. 

“Enhancement of band strength f Sa{X)d\~^ relative to average of lab values (see Table 1 of LDOl). 

‘^Parameters recommended by LDOl for diffuse regions in the Milky Way. 

“Parameters giving good agreement with position and FWHM of 5-15/rm features for the SMC Bl^l molecular cloud 
(see §4.2). 

^Band strengths giving good agreement with spectrum of the SMC Bl^l molecular cloud (see §4.5). 


The PAH absorption cross sections Cabs (o; A) and the temperature probability distribution 
functions dP/dT depend upon the number of C atoms and the H/C ratio in the molecule. As a 
starting point, we take H/C to be 


H/C 


' 0.5, Ac < 25, 

< 0.5/v^Ac/25, 25 < Ac < 100, 
0.25, Ac > 100, 


( 11 ) 


which is appropriate for compact, symmetric PAHs for Nq < 10^, and clusters of such PAHs for 
Ac > 10^. Draine & Li (2001) found that the prescription (11) appears to be consistent with the 
3-15^m emission spectrum seen in reflection nebulae and photodissociation regions in the Milky 
Way. 


For a given u-e and a given radiation intensity A, we estimate the PAH ionization fractions 
0ion(A, a) in SMC Bl^l using the rates for photoelectric emission and electron capture recently 
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discussed by Weingartner &: Draine (2001b). Figure 3 shows the ionization fractions calculated for 
a range of radiation intensities and PAH sizes. At low U/ue-, PAHs are more likely to be negatively 
charged (e.g., (/>ion+ < l^ion- U/ne ^ 10 cm^; see Figure 3a,b) since the charging of a neutral 
PAH is dominated by electron capture. 


20 50 102 


Nc 

102 


10 “ 



Fig. 3.— The fraction of positively charged (a), negatively charged (b), and non-neutral (c) PAHs as a function 
of size and radiation intensity for Tie = 0.1 cm“® and Tgas = 15 K. The upper axis label gives the number of carbon 
atoms Nc = 0.468(a/ A)®. 

For each of our two trial Ue values, we adjust ?7miin tAnaxj /3, o-o, and be to minimize x^- In 
Figure 4a,b we plot the best-fitting model spectrum for our low-ng model, for which t/min = 

Uma,x = 10^'^, /? = 2, oq = 3A, a = 0.4, be = 1.5ppm, and 4.2 x 10^^ cm“^ (parameters 

are summarized in Table 1). Our estimate for the mean radiation intensity {U) ~ 2.2 is consistent 

with 17 ~ 3 (with an uncertainty of a factor 2) estimated by Reach et al. (2000). We estimate that 

the dust in SMC provides an extinction Ay ~ 2.1 mag. 

Although the model closely reproduces the 12, 25, 60, 100, 140, and 240/xm IRAS/DIRBE 
photometry (the 2.2, 3.5, and 4.9//m DIRBE bands are well represented by 6500 K black-body stellar 
radiation), the model is not fully successful in reproducing the 6-13/im PAH emission features, 
resulting in the relatively large x^/d.o.f. ~ 5.3. In particular, it is a bit too weak for the 11.3/im 
C-H out-of-plane bending mode while too strong for the 7.7^m C-C stretching mode. 

In view of the fact that the 11.3/im feature is much stronger in neutral PAHs than in charged 
PAHs (Allamandola, Hudgins, & Sandford 1999), the 11.3/im feature deficit could be alleviated by 
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increasing ng so as to reduce the PAH ionization fraction (see Figure 3). Indeed, it is seen in Figure 
4d that our high Ue model provides an improved fit, with x^/d.o.f. ~ 3.1. The unphysical case of 
pure neutral PAHs leads to a better fit (see Figure 4f; y^/d.o.f. 1.9), but is still deficient in the 

11.3^m band. We discuss this further in §4.3. 

For both low-rig and high-ng models, the number of C atoms (relative to H) in PAHs required 
to account for the SMC mid-IR spectrum is he ~ 1.5ppm. The total C and Si abundances 

locked up in dust ([C/H]dust ~ 24.2 ppm; [Si/H](just ~ 4.1 ppm) are consistent with the overall SMC 
elemental abundance (C/Hfs 54 ppm, Si/H^i 11 ppm, Mg/Hss 9.6 ppm, Fe/H^i 6.9 ppm; Russell & 
Dopita 1992). 

Finally, we remark that models with a single value of U are unable to provide a satisfactory 
fit to the SMC V>\^1 60^m flux unless we add a log-normal PAH component peaking at ~ 40 — 
50A containing a substantial mass (e.g. C/H > 8 ppm for U = 3). 



Fig. 4.— Comparison of the model (solid line) to IRAS photometry (squares), DIRBE photometry (diamonds), 
and 5-15/im ISOCAM spectroscopy (vertical bars in panels b, d, and f) of the quiescent molecular cloud SMC Bl^l. 
Triangles show the model spectrum convolved with the IRAS (12, 25, 60, lOOrtm) and DIRBE (140, 240rim) filters. 
Stellar radiation (dot-dashed line) is approximated by a 6500 K black-body. Upper panels (a,b): the best-fit low-rie 
dust model (ao = SA, cr = 0.4, be = 1.5 ppm; rie = 0.1 cm”®; x^/d.o.f. « 5.3) illuminated by a range of radiation 
intensities dNu/dU oc U~^, 10“° ® < U < 10®'®, with ss 4.2 x 10^^ cm”®. Middle panels (c,d): same as (a,b) but 
for the high-rie model [ue = 1.0 cm”®; x®/d.o.f. ~ 3.1). Lower panels (e,f): same as (a,b) but for a pure neutral-VKH 
dust model (x®/d.o.f. Ri 1.9). 
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4.3. Models with Increased PAH Hydrogenation 

It is seen in Figure 4 that even when entirely neutral our PAH model is still too weak in the 
11.3//m band. Reach et ah (2000) suggested that the strong emission in the 11.3/im C-H bending 
mode might be due to increased hydrogenation of the SMC PAHs, which might perhaps be related 
to the fact that the gas phase C abundance (relative to H) is a factor of 10 lower than in the 
Milky Way (Dufour 1984; Russell & Dopita 1992). To test this hypothesis, we have considered a 
PAH model where we assume the same intrinsic band strengths for C-C and C-H stretching and 
bending modes, but with a higher H/C ratio (which would be appropriate for PAHs with more 
open structures); 


H/C 


' 0.5, Nc < 100, 

< 0.5/y^Ac/100, 100 <Nc< 400, 

0.25, Nc > 400. 


( 12 ) 


We find that the best-fitting spectra (see Figure 5) do lead to improvements over those calculated 
from PAHs with a lower degree of hydrogenation (see Eq.[ll]). 


However, this high H/C model still fails to attain the observed 11.3/im/7.7/im band ratio even 
taking all PAHs to be neutral (see Figure 5c). An increase in oq and/or a (implying a larger c/ion 
[see Figure 3] and a smaller H/C [see Eq.l2]) does not help since this would lead to a decrease at 
the 11.3;um band. While a slight decrease in oq and/or a (implying a smaller c/ion [see Eigure 3] 
and a higher H/C [see Eq.l2]) would result in an increase at the 11.3/im band, this results in a 
decrease at the 6.2, and 7.7^m bands and thus does not work either since this model is already 
deficient in the 6.2 and 7.7^m bands (see Figure 5f). 


In comparison with normal PAHs, super-hydrogenated PAHs would produce a stronger 3.3^m 
C-H stretching feature, a broad plateau below the 3.4 and 3.5^m features (Bernstein, Sandford, & 
Allamandola 1996; Pauzat & Ellinger 2001), and a stronger 6.2^m C-C stretching feature (Beegle, 
Wdowiak, & Harrison 2001). Future high-resolution observations of SMC particularly at 

the 3^m region, will provide us with more information about the PAH hydrogenation. 


Since Milky Way PAHs with >25 C atoms are expected to be already essentially fully hy¬ 
drogenated (Tielens et al. 1987; Allamandola, Tielens, & Barker 1989; Allain, Leach, & Sedlmayr 
1996), an increase in the H/C ratio in the SMC would require that the PAHs have more open 
structures than in the Milky Way. There is no apparent reason to expect this difference in PAH 
structure between the Milky Way and SMC. Based on an analysis of the 10-15/um PAH out-of-plane 
C-H bending features, Vermeij et al. (2002) argue that the PAH emission in the SMC Bl^l cloud is 
dominated by compact PAHs, i.e., relatively low hydrogenation. In conclusion, it appears unlikely 
that the large 11.3;um/7.7/um band ratio in SMC Bl^l is attributable to enhanced hydrogenation. 



14 



Fig. 5.— Same as Figure 4 but for PAHs with a higher degree of hydrogenation (see Eq.[12]). 

4.4. Reddened Starlight? 

In previous sections the interstellar radiation field is taken to have the same spectrum as 
the MMP solar neighbourhood field. Although this assumption seems reasonable for the overall 
SMC where the extinction is small {Ay ~ 0.3mag; see §3), the starlight incident on the SMC Bl^l 
molecular cloud surface will be subject to considerable reddening inside the cloud {Ay ^ 2 mag; see 
§4). The reddened starlight will affect the PAH emission bands in two ways: (1) PAH molecules will 
be excited to lower temperatures; (2) the probability of being (positively) charged (via photoelectron 
ejection) will be smaller. One therefore would expect an increase in the 11.3/xm/7.7/xm band ratio. 
Could reddening of the starlight explain the observed 11.3/im/7.7/im band ratio in SMC 

To investigate this, we consider a trial reddened radiation field with the spectrum shown in 
Figure 6a, with a factor of ~3 reduction in the relative intensity at O.l^m relative to 0.5/rm. We have 
calculated the temperature probability distribution functions for PAHs excited by this radiation 
field with (7 = 1. The mid-IR emission spectrum is entirely due to single-photon heating, so the 
emission per PAH simply scales as U ; the spectrum in the lower panel of Figure 6 is calculated 
for U = 2.2. It is seen in Figure 6b that even for pure neutral PAHs, reddened starlight results in 
an emission spectrum which is very similar to the pure neutral PAH spectrum in the unreddened 
radiation field in Figure 4f, and is unable to reproduce the observed 11.3/im/7.7/im band ratio. We 
conclude that excitation by reddened starlight cannot account for the 11.3^m/7.7^m band ratio 
observed toward SMC Bl#l. 
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Fig. 6 .— Upper panel (a): trial reddened ISRF spectrum (solid line) with energy density equal to the MMP ISRF 
[U = 1). Also plotted is the MMP ISRF {U = 1; dashed line). Lower panel (b): the best-fit pure neutral-PAH model 
(oq = 2.6 A, cr = 0.4, be = 1.5 ppm, and A^h ~ 4.5 x 10^^ cm“^ [for {U) ~ 2.2); y^/d.o.f. ~ 1.9). 

4.5. IR Band Strengths: Milky-Way PAHs vs. SMC PAHs 

Individual PAHs differ considerably in the relative strengths of different vibration bands. LDOl 
discussed laboratory measurements and theoretical calculations of band strengths for selected PAH 
molecules and ions, but argued that spectroscopy of various objects in the Milky Way was best 
reproduced if the 6.2, 7.7 and 8.6/im band strengths were increased by factors of Eq ,2 = 3, FCy.y = 
2, and Flg.g = 2 respectively, relative to the average laboratory and theoretical values. These 
enhancement factors were used for the models shown in Figures 4-6. 

It is reasonable to ask whether the observed spectrum indicates that the band strengths for 
the PAHs in SMC Bl^I differ from the band strengths which characterize PAHs in the Milky Way. 
If so, what band strengths give the best fit to the SMC BI^I spectrum? 

In Figure 7 we show the emission calculated from a model with £’ 5.2 = 1-5, £ 7.7 = £§.6 = I- 
For He = 1.0 cm“^ this model (see parameters in Table 1) is in very good agreement with the 
observations, with y^/d.o.f. = 1.7. The agreement is considerably better than for the models 
considered in Figures 4-6. Note that these adjusted band strengths are in fact closer to the average 
laboratory values than are the “Milky Way” band strengths. 

Alternatively, the apparent need for a change in intrinsic band strength might be due to the 
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oversimplification of the “astronomical” PAH model in which we do not distinguish anions (PAH”) 
and multiply-charged cations (PAH”'"*') from singly-charged cations (PAH"*"). Bauschlicher &: Bakes 
(2000) predict that the 6.2, 7.7//m C-C stretching modes and the 8.6//m in-plane C-H bending mode 
are even more enhanced for PAH”"*" than PAH"*". PAH” anions have band strengths intermediate 
between those of neutrals (strong 3.3^m C-H stretching and 11.3, 11.9, 12.7/rm out-of-plane C-H 
bending modes) and PAH"*" cations (strong 6.2, C-C stretching and 8.6/xm C-H in-plane 

bending mode). Unfortunately, no experimental measurements have yet been made for multiply- 
ionized PAHs. Available evidence indicates that the IR properties of PAH” anions closely resemble 
those of PAH"*" cations except for the very strong 3.3//m C-H stretch enhancement in the anion 
(Szczepanski, Wehlburg, &: Vala 1995; Hudgins et al. 2000). 



Fig. 7.— Same as Figure 4 but for PAHs with “SMC” band strengths (see text) and a size distribution of ao=2.8A 
and a—0.4. We favor this model with rig ~ 1.0 cm”® for SMC Bl^l. 


4.6. UV Extinction 

Laboratory measurements show that PAHs have an absorption peak due to the tt-tt* electronic 
transition in the vicinity of 2000-2500A, and in the grain model considered here (LDOl) it is assumed 
that the interstellar PAH mixture contributes absorption with the wavelength-dependence of the 
2175 A feature observed on sightlines in the Milky Way and the Large Magellanic Cloud (LMC). 
The PAH abundance required to reproduce the observed IR emission features in the Milky Way 
(C/H~ 60 ppm) is such that we posit that PAHs in fact contribute most of the observed 2175A 
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feature. The predicted extinction excess at 2175A is (see LDOl) 

^^2175 A ~ {bcNn/W^'^ C cm"^) mag . (13) 

Our best-fitting model has A^h = 4.2 x 10^^ cm“^ and be = 1.5 ppm (see §4). Therefore we predict 
AT » ~ 0.5 mag for sightlines through SMC Bl^l. We note that the AT « value derived 
here is a lower limit since 15 < a < 150 A carbonaceous grains may also contribute to the 2175A 
hump. 

This prediction appears to conflict with the fact that most SMC extinction curves have no 
detectable 2175A hump (Prevot et al. 1984). However, as pointed out by Reach et al. (2000), 
there may be regional variations in the SMC extinction curve. We note that the sight lines which 
show no 2175 A hump all pass through the SMC Bar regions of active star formation (Prevot et 
al. 1984; Rodrigues et al. 1997; Gordon & Clayton 1998). There is at least one line of sight (Sk 
143=AvZ 456) with an extinction curve with a strong 2175A hump (Lequeux et al. 1982; Prevot 
et al. 1984; Bouchet et al. 1985; Thompson et al. 1988; Rodrigues et al. 1997; Gordon & Clayton 
1998). This sightline passes through the SMC wing, a region with much weaker star formation 
(Gordon & Clayton 1998). It is possible that the SMC Bl^l dust, just like the dust toward Sk 
143, has been exposed to a less harsh environment than that in the SMC Bar regions where the 
UV extinction has been measured, and where the radiation and shocks associated with massive 
stars might have destroyed PAHs. Finally, we also note that both the 2175A hump (Fitzpatrick 
1985) and PAH emission features (Sturm et al. 2000) are present in the 30 Dor region in the LMC, 
another low-metallicity environment (~l/4 of that in the Milky Way; Dufour 1984). 

5. Conclusions 

In contrast to the Milky Way Galaxy, the SMG is characterized by a lower metallicity, a lower 
dust-to-gas ratio, and an absence of the 2175A hump on most sightlines with measured UV extinc¬ 
tion (see §1). Therefore, the SMG provides an ideal laboratory to study dust in an environment 
quite different from the Milky Way. We have modelled the IR spectral energy distribution for the 
SMC as a whole (averaged over a 6.25 deg^ area) (§3). We have also modelled the IR emission 
spectrum of the quiescent molecular cloud SMC Bl^l (§4). Our principal results are: 

1. The dust IR emission from a 6.25 deg^ region including the SMC Bar and Eastern Wing 

is well reproduced by a model with dNn/dU oc 10“^'° < U < 10^'^®, and ^ 

5.4 X 10^^ cm“^ using the SMC dust model of Weingartner & Draine (2001a). This dust has a 
very low PAH abundance, with C/H < 0.2ppm (< 0.4% of the SMC interstellar C abundance) 
incorporated into PAHs. The dust-weighted mean radiation intensity is {U) ~ 2, but the 
starlight intensity must have a broad distribution ranging from Umin ^ 0.3 to Umax ^ 50 (§3). 

2. The IR emission from the molecular cloud SMC Bl^l can be reproduced by a dust model 
consisting of silicates, graphite, and PAHs with ~ 4.2 x 10^^ cm“^ illuminated by starlight 
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with a power law distribution dNu/dU oc. U ^,10 < U < 10^'^ (§4). The required PAH 

abundance in SMC Bl^l is 1.5 ppm, or ~3% of the SMC interstellar C abundance. 

3. The observed 11.3^m/7.7/um and 6.2/rm/7.7/xm band ratios for SMC Bl^^l fall outside the 
“allowed” region predicted for the Milky Way PAH model (see Figures 16, 17 of Draine & Li 
2001). In §4.2, §4.3, and §4.4, we examine the possibility that this could be due to (1) an 
enhanced neutral fraction, (2) an enhanced H/C ratio for the PAHs in SMC Bl^l, or (3) 
reddening of the radiation field to which the PAHs in SMC Bl^)!! are exposed, but none of 
these scenarios successfully accounts for the observed 5-15/um spectrum. We instead conclude 
that the intrinsic IR band strengths for the PAH mixture in SMC Bl^l differ from the band 
strengths adopted for Milky Way PAH mixtures (§4.5). This indicates that the PAH mixture 
in SMC Bl^l differs in detail from the PAH mixtures present in a number of Milky Way 
regions. The band strengths suggested by the SMC Bl^^tl spectrum are close to “laboratory” 
band strengths. 

4. The PAH abundance in SMC BItt^:! exceeds that in the SMC Bar by a factor > 8 (C/H=1.5 
ppm vs < 0.2ppm). While the average SMC extinction curve has no detectable 2175A hump, 
we predict a 2175A extinction excess of > 0.5mag for the dust in SMC Bl^l (§4.6). 

Two types of measurements would be of great value to test our dust model; (1) UV extinction 
measurements for sightlines through or near SMC Bl^l to test for the predicted 2175A feature; 
(2) 3 — 15^m spectrophotometry of other regions in the SMC where UV extinction measurements 
place a strong upper limit on the 2175A feature - the present model predicts that PAH emission 
should not be seen from such regions. 


We thank W.T. Reach for providing us with the ISO spectrum of SMC and for helpful com¬ 
ments; D.P. Finkbeiner for help in obtaining the DIRBE fluxes for the SMC BI^t^I molecular cloud; 

S. Stanimirovic for obtaining the SMC HI mass; G. Helou for valuable comments; and R.H. Lupton 
for the availability of the SM plotting package. This research was supported in part by NASA grant 
NAG5-10811 and NSF grant AST-9988126. 
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